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Abstract 

The collective behavior in a variant of Schelling's segregation model is characterized with 
methods borrowed from statistical physics, in a context where their relevance was not conspic- 
uous. A measure of segregation based on cluster geometry is defined and several quantities 
analogous to those used to describe physical lattice models at equilibrium are introduced. This 
physical approach allows to distinguish quantitatively several regimes and to characterize the 
transitions between them, leading to the building of a phase diagram. Some of the transitions 
evoke empirical sudden ethnic turnovers. We also establish links with 'spin-1' models in physics. 
Our approach provides generic tools to analyze the dynamics of other socio-economic systems. 

Keywords: segregation - Schelling model - phase diagram - discontinuous transition - 
Blume-Capel model - Blume-Emery-Griffiths model 

Introduction 

In the course of his study of the segregation effects observed in many social situations, Thomas 
Schelling introduced in the 1970's [T] [2] a model that has attracted a lot of attention ever since, 
to the point that it may now be considered an archetype in the social sciences. The success of 
the Schelling model is due to several factors: It was one of the first models of a complex system 
to show emergent behavior due to interactions among agents; it is very simple to describe, yet its 
main outcome - that strong segregation effects can arise from rather weak individual preferences 
- came as a surprise and proved robust with respect to various more realistic refinements; as a 
consequence, it has possibly far-reaching implications for social and economic policies aiming at 
fighting urban segregation, considered a major issue in many countries (for recent discussions of the 
social relevance of the model, see [5] and [1]). 

An important recent development is the realization that there exists a striking kinship between 
this model and various physical models used to describe surface tension phenomena [5 or phase 
transitions and clustering effects, such as the Ising model 00IE]- This connection is not a rigorous 
correspondence, still it is more than a mere analogy. It gives novel insight into the behavior of 
the Schelling model, and it suggests more generally that socio-economic models may be fruitfully 
attacked drawing from the toolbox of statistical physics [Tj |9] ■ Indeed, physicists developed during 
the last decades powerful methods for situations where obtaining analytical results seems out of 
reach. These rely on the quantitative analysis of computer simulation results, guided by some 
general principles. They are well suited to complex systems such as those encountered in the social 
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sciences and should in particular prove powerful in conjunction with agent-based modeling 1 , a 
growingly popular approach |10j . 

A key physicist strategy is to characterize the system under study by the phase diagram which 
gives, in the space of control parameters, the boundaries separating domains of different qualitative 
behaviors. Each type of behavior is qualified by the order of magnitude of a small set of macroscopic 
quantities, the so-called "order parameters" . The main difficulties are to correctly identify the 
relevant set of order parameters and the associated qualitative behaviors of the system, and to locate 
and characterize the boundaries, on which "phase transitions" occur in a smooth or discontinuous 
way. In the present paper, we show how several of the methods evoked above can be adapted and 
applied to the building of the phase diagram of social dynamics models, taking as paradigm the 
Schelling segregation model. More precisely, we illustrate our approach on a particular variant of 
the Schelling model, where the basic variables are the tolerance (to be defined precisely below) 
and the density of vacant sites. We introduce two order parameters which provide a relevant 
measure of global segregation, a surrogate of the energy, and analogues of the susceptibility and 
the specific heat. We also introduce a real-space renormalization method suited to situations with 
a high density of vacancies. Our analysis shows the existence in the phase diagram of the model 
of sharp transitions, where relevant quantities have singularities, and we discuss the nature of 
these transitions. They separate different types of behavior - segregated, mixed, or "frozen" -, in 
agreement with qualitative observations based on pictures of simulated systems [5] . In addition we 
make contact with the Blume-Emery-Griffiths jTT] and Blume-Capel [T5J[T3] models, which are spin- 
1 models much studied in relation with binary mixtures containing mobile vacancies and which show 
a richer behavior than the simple Ising model, with both discontinuous and continuous transitions. 
One should insist that most variants of the Schelling segregation model are of kinetic nature: their 
dynamics cannot be described as the relaxation to an equilibrium characterized by some energy 
function, except for specific variants (for such exceptions, see e.g. |14l 121] and the discussion below 
on the links with spin models). Nevertherless, the tools and quantities we introduce by analogy with 
equilibrium statistical mechanics appear to be quite efficient to characterize the model behaviors. 
They are sufficiently general and relatively simple to be adaptable to a large variety of social and 
economic models, as long as these involve interacting agents living in a discrete space or more 
generally on a social network. 

1 Model and qualitative analysis 

In Schelling's original model [I] agents of two possible colors are located on the sites of a chessboard. 
Each color corresponds to members of one of two homogeneous groups which differ for example by 
their race, their wealth, etc... A fraction of the sites are blanks, the agents of both colors may move 
to these vacancies. The neighborhood of an agent comprises the eight nearest and next-nearest sites 
(Moore neighborhood). If less than | of an agent's neighbors belong to his group, he is discontent 
- in economic terms his utility is 0; otherwise he is satisfied - his utility is 1. Starting from random 
initial configurations Schelling displaced discontent agents onto the closest satisfactory vacant sites, 
if possible. He observed that the system always reached a segregated state, where large clusters 
of same-color agents were formed. The crucial point is that segregation appears as an emergent 
phenomenon, in the sense that the collective effect is much stronger than what would be naively 
expected, as individual agents are happy to live in a mixed neighborhood. This phenomenon proves 
robust: a similar outcome, with some caveats, is found in variants of the model, even when the 
utility function is non-monotonous with the fraction of similar neighbors [51 115] . 
The model we consider is a variant of the original Schelling model: The agents are satisfied with 
their neighborhood if it is constituted of a number of unlike agents lower than (or equal to) a 
fixed proportion T of all the agents in the neighborhood. The parameter T is called the tolerance 
[I] . Since a higher tolerance allows for a larger number of configurations of satisfied agents (higher 
entropy), this parameter may be thought of as a temperature-like variable. We will be guided by 
this qualitative correspondence - although in a different way than in [6 , where a direct analogy is 

lr Thc 'hand-made simulations' done by T. Schelling by moving pawns on a chessboard can be considered as the 
first agent-based simulations ever done in social science. 
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made with the Ising model. The other control parameter of the model is the vacancy density p. The 
randomly chosen agents move one by one to any vacancy which has a satisfying neighborhood - this 
is equivalent to long-range diffusion in physical terms. If no vacancy fits for some agents the latter 
respectively move back to their initial position. This dynamics is repeated until configurations are 
reached where the number of satisfied agents is almost stable. Note that in the present variant 
satisfied agents can also move, not just discontent ones. That rule introduces some noise in the 
dynamics and is useful to avoid a particularity of the original Schelling model noted in [6l [14] , 
namely that the system may end up in states where the clusters are large but finite, so that strictly 
speaking no large-scale segregation occurs. We will see later how the intensity of this noise is actually 
correlated with the tolerance level T. Note also that the global utility may decrease at times during 
the process, as the gain for the moved agent can be less than the net loss for his old plus new 
neighbors [2 ]. Let us emphasize that only a finite number of values of T are meaningful, namely 
i, i..., |, |. They correspond to the maximal number of tolerated different neighbors divided by 
the actual number of occupied sites in the neighborhood. Any other value for T is thus equivalent 
to the closest inferior meaningful value. 

1.1 Numerical Simulations. 

Our simulations were performed on a L x L lattice (L = 50, unless otherwise specified) with free 
boundary conditions. An initial configuration was randomly generated such that the vacancies and 
the two types of agents were fully mixed. Then the evolution followed the rules described above. 
In the simulations one time step corresponds to one attempted move per occupied site on average, 
the usual definition of a Monte-Carlo step. 




St=0 St=10 St=100 St=400 



Figure 1: Evolution of the configuration for a vacancy concentration p = 5% and a tolerance T = 0.5 
with a network size L = 100. St stands for the number of time steps. The red and blue pixels 
correspond to the two types of agents, the white pixels to the vacancies. The system evolves from 
a random configuration - where the vacancies and the two types of agents are intimately mixed - 
to a completely segregated configuration. After just 10 steps there exist two percolating clusters, 
one of each color, which are very convoluted, fractal-like. 



Fig. [T] shows the time evolution of a typical configuration, for a vacancy density of 5% and a 
tolerance T = 0.5, which means that agents accept at most half of their neighbors to be different 
from themselves. One observes the rapid formation of large clusters. After 10 time steps, the 
proportion of satisfied agents is already very close to 1 (Appendix, Fig. A.2|. It increases slowly 
thereafter, but the structure of the clusters keeps evolving. They become more and more compact 
and well separated spatially, their surface gets less corrugated, in a process strongly reminiscent 
of coarsening effects in alloys [TB] • We now consider more general values of the tolerance and the 
vacancy density. Fig. [2] shows configurations obtained after letting the system evolve with the 
dynamics previously described until it reaches equilibrium. What is meant by equilibrium here may 
correspond to two different situations: (i) The system does not evolve at all anymore (fixed point); 
(ii) The systems reaches some stationary state: the fluctuations of the studied parameters remain 
weak during a large number of time steps. In the following all averaged quantities are measured 
during 30000 steps after equilibrium is reached. 

At small and moderate values of p, one observes that : 

- For low values of the tolerance the system stays in a mixed state, no large one-color clusters 
are formed although this would be more satisfactory for the agents. Actually, whatever the initial 
configuration, the system remains close to the state in which it was created: this is a dynamically 
frozen state. 
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Figure 2: Configurations obtained at large times for selected values of p and T. 



- When T increases, at fixed p, a drastic qualitative change occurs at an intermediate value of 
T: The system separates into two homogeneous regions of different colors, segregation occurs. This 
behavior subsists for an interval of T which depends on p. 

- For T larger than a value weakly depending on p, the final configuration is again mixed. 
For large values of p one observes a smooth transition as T increases from a segregated state to 
a mixed one. We will now give a more quantitative description and characterize the transitions 
between these different behaviors. 



2 Quantitative analysis: order parameters 
2.1 Main order parameter: A measure of segregation 

Though the presence of segregation in a system can be visually assessed, a quantitative way to 
measure it is necessary, in particular to discuss the nature of the transitions between different 
states. Different possible measures have been suggested, by Schelling himself and later by various 
authors [3lll4[IT7]. which capture various aspects of the phenomenon. Here we introduce a measure 
linked to the definition of segregation as the grouping of agents of the same type and the exclusion 
of the other type in a given area. To that effect we consider that two agents belong to the same 
cluster if they are nearest neighbors (the Moore neighborhood cannot be used here, since Moore 
clusters of different colors may overlap and be large without segregation occurring) . The (mass) size 
of a cluster c, i.e the number of agents it contains, will be called n c . Taking a hint from percolation 
theory where it plays a central role [15] . we introduce the weighted average S of the size of the 
clusters in one configuration 

S =) j n c p c (1) 
M 

where p c = is the weight of cluster c, N to t — L 2 (\ — p) being the total number of agents. The 
maximal size of a cluster is N tot /2, so the normalized weighted cluster size s is given by 

s = U{hpj s= ( L2 (1 2 _ p))2 I> 2 ( 2 ) 

The sample average of s after reaching equilibrium will be called the segregation coefficient (s) . 
Its value for complete segregation (i.e., only two clusters survive) is 1 and it vanishes if the size 
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of the clusters remains finite when the system dimension L tends to infinity. It may therefore 
play the role of an order parameter to identify a segregation transition. The variation of the 
segregation coefficient (s) with respect to the tolerance is illustrated in Fig. [3] for different values of 
the vacancy density. The calculations were done using the Hoshen-Kopelman algorithm to labelize 
the clusters [TH]- For each density of vacancies there exist two critical values of the tolerance. At 




0.2 0.4 0.6 0.8 
Tolerance T 

Figure 3: Segregation coefficient (average of s defined in Eq. ^) for several values of the vacancy 
density p. The lines linking the points are guides to the eye. 

the first one Tt (/ for frozen), (s) jumps from a very low value to about 1. This signals an abrupt 
change from a mixed configuration to one with only two clusters (one for each type of agent). A 
second jump, in the reverse direction, occurs for a larger tolerance T c . This second value depends 
slightly on p. unlike Tf. The higher the density of vacancies, the smaller the value of T at which 
the segregation phenomenon appears and the broader the interval of T for which it exists. For a 
vacancy density above about 20% the segregation coefficient departs from 1 (Fig. |3|, even if the 
agents are visually segregated (Fig. |2.1| . This is due to the definition of a cluster based only on the 
four nearest neighbors. In some regions, even if there are agents of one color only the presence of 
many vacancies may lead to group them into distinct clusters and to miss the existence of a "dilute 
segregation" situation. In order to identify clusters at a larger scale in such cases we introduce 
a real-space renormalization |20j procedure. An example of the renormalization process for a site 
and its neighborhood is illustrated on the left of Fig |2.1| A renormalized configuration is shown 
on the right part of Fig|2.1| (see the appendix for details on the procedure) . The renormalization 




Figure 4: Example of renormalization. The renormalization is performed on a configuration corre- 
sponding to T = g and p — 50%. 

has a strong effect for high values of the density of vacancies. For p = 50% the configurations are 
visually segregated for a range of small and medium values of the tolerance but the raw segregation 
coefficient is very small (Fig. [5] left), whereas after renormalization it is very close to 1 (Fig. [5j 
right). We can conclude that, for values of the vacancy density strictly less than 50%, there is a 
discontinuous transition from the segregated to the mixed state (as shown on Fig. H for p up to 
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Figure 5: Segregation coefficient as a function of the tolerance for p — 50%: Raw data (left) and 
results for the renormalized systems (right). These suggest that a continuous transition takes place 
near T = 0.6. The error bars were obtained by computing the variance a s , a 2 = (s) 2 — (s 2 ) (same 
for Figj3j. 



26%). At p = 50% the transition becomes continuous - as shown on Fig (5] the order parameter 
has a smoother variation -, and there is no longer any clear transition for vacancy densities above 
~ 56% (not shown). Moreover, the configurations observed at large p suggest the existence of a 
diluted phase of segregation at high vacancy densities: the two types of agents are not mixed but 
there are domains with many small clusters of a same color in a sea of vacancies. As for the frozen 
state observed at low T for small p, it disappears at some medium value of p (compare Fig. [3] with 
Fig. [5j right). 

2.2 A second order parameter: densities of unwanted locations. 

The segregation coefficient defined above does not allow to distinguish between the mixed state at 
low tolerance and the one at high tolerance (Figj3j). However, even if the final configuration is a 
mixed configuration for both situations, the nature of the two states is not the same. An additional 
parameter is thus necessary to analyze the results. We introduce the density p r of empty places 
where the red agents do not want to move (symmetrically p~b for blue agents) . The plot of p r versus 




Figure 6: Density p~ r of vacancies where the red agents would be unsatisfied, for several vacancy 
densities. The results for the blue agents are similar. Indeed, they play symmetrical roles in the 
present model. 
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the tolerance (Figj6J) shows that this quantity undergoes two jumps for each vacancy density p. 
Until the tolerance reaches X 1 / no empty space is considered attractive by these agents. Between Tf 
and T c half of the vacant sites are satisfactory for one type of agents. For T larger than T c almost 
all empty spaces are acceptable. A representative configuration for a situation where all moves are 
permitted is an unordered one, this explains the mixed situation observed at high tolerances. The 
quantity p r thus allows to characterize the three regimes and to discriminate between the low and 
high tolerance mixed situations. 

Thanks to the quantities defined previously we have identified the regions of segregation. Let 
us summarize the main results obtained so far. The final configurations are reminiscent of those 
encountered in spin lattice models with paramagnetic (= mixed state) and ferromagnetic (= segre- 
gated state) phases. The segregation coefficient and the density of unilaterally unwanted locations 
play the role of order parameters. They show discontinuous jumps at some particular values of the 
tolerance T, indicating the existence of sharp phase transitions. For not too large values of the 
vacancy density there are two such transitions: At low tolerance, the system goes from a frozen 
state due to the dynamics that does not allow any movement, to a segregated state. At higher 
tolerance the system becomes mixed again. 

3 Analogues of thermodynamic quantities 

3.1 Contact with Spin-1 models and analogue of the energy. 

We now introduce several other useful quantities by analogy with thermodynamic properties studied 
in statistical physics. To do so, we first exhibit a link with spin-1 models. 

One can show that (see the appendix and [21]), had we forbidden the displacement of satisfied 
agents (as studied in [21 017]), the dynamics would have a Lyapunov function, that is a quantity 
which decreases with time, driving the system towards a fixed point. This function can be written 
under the form 



where K = 2T — 1 and the CjS are 'spin-1' variables taking the value if the location i is not 
occupied and 1 (resp. —1) if this location is occupied by a red (resp. blue) agent; the sums are 
performed on the nearest and next nearest neighbors. This function ^ is identical to the energy 
of the Blume-Emery- Griffiths model [11 under the constraint that the number of sites of each type 
(0,±1) is kept fixed. This spin-1 model, and the Blume Capel model [T^llIS] corresponding to the 
particular case K = 0, have been used in particular to modelize binary mixtures and alloys in the 
presence of vacancies. A more detailed analysis of the link between the Schelling model and the 
Blume-Emery- Griffiths model will be published elsewhere [21 . In the particular variant considered 
here, however, the dynamical rules do not lead to the minimization of such a global energy. Yet, it is 
clearly potentially interesting to consider this quantity Es as a surrogate of the energy. Compared 
to the dynamics having Es as a Lyapunov function, the moves of satisfied agents introduce a source 
of noise which has some similarity with a thermal noise. Its amplitude may be measured by the 
fraction of agents who are satisfied. When starting from a random initial configuration this is higher 
at higher tolerance, hence the tolerance value is an indirect measure of this noise level. This gives 
another motivation, different from the one already evoked, for taking the analogy between T and a 
temperature as a guideline for the analysis, as done in what follows. We find that the average of the 
second part of the energy Es essentially consists of a term linear with K (see appendix) , so that 
the transitions are more easily located by only plotting the average of the first term of Es, Fig. [7] 
corresponding to the Blume-Capel part of the energy, Ebc — ~ "Ylu j) c i c j- It confirms the existence 
of the two transitions previously evoked: At low tolerance its decrease occurs at the transition from 
the frozen state to the segregated one, whereas the increase observed at high tolerances corresponds 
to the transition to the mixed state. Such abrupt variations are characteristic of a discontinuous 
- in thermodynamic language "first order" - transition. Moreover we note from Fig. [7] that as the 
vacancy density increases the energy varies less abruptly: This signals a change in the nature of the 
transition, from discontinuous to continuous ('second order'), as will be discussed below. 
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Figure 7: Variation of the mean of the Blumc Capel energy Esc for different values of the vacancy 
density. The data are normalized by 4L 2 (1 — p). 



3.2 Analogue of the specific heat 

Since the energy analogue proved fruitful, an analogue of the specific heat may also be expected 
to give useful information. However, here we cannot make use of the thermodynamic definition 
C s = dEg/dT since only a finite number of values of the tolerance have physical meaning. In order 
to work with a well-defined quantity at fixed tolerance T we remark that the specific heat is related 
to the energy fluctuations at equilibrium via the so-called fluctuation-dissipation theorem. In the 
present context the relevant formulation of this important theorem is: 

C=<S>^ (4) 

where the notation ( ) means an average over the configurations taken by the system after reaching 
cquibrium, and T is the tolerance. We will call C s the fluctuation coefficient. It plays a role analo- 




Figurc 8: Variation with the vacancy density of the fluctuation coefficient C s . 



gous to the volatility index measuring price fluctations in financial markets. 

According to Fig. [8] the fluctuation coefficient has a well-marked peak at the "segregated-mixed" 
transition. This peak flattens out and tends to disappear as the vacancy density increases, confirm- 
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ing the disappearance of the transition from segregated to mixed states. On the other hand this 
fluctuation coefficient, null at very low tolerance, has a slight increase at the "frozen-segregated" 
transition. The jump in the energy Eg does not show up in the fluctuation coefficient C s , indicat- 
ing that the fluctuation-dissipation theorem is stongly violated, as typically observed in dynamic 
transitions in glassy or kinetically constrained systems |23j . 

3.3 Analogue of the susceptibility 

In the same vein one would like to define the susceptibility, a quantity which is to the order- 
parameter what the specific heat is to the energy. Here, having proposed the segregation coefficient 
as the main order parameter, we define the analogue of the susceptibility through the fluctuation- 
dissipation relation involving the order parameter fluctuations 



Xs 



<, 2 



T 



(5) 



where s is defined by Eq. ^ . As can be seen on Fig. [9] this susceptibility presents a peak at the 
second transition which is more conspicuous for high vacancy densities than the one for C s . This 
peak does not exist anymore for high vacancy densities in agreement with the disappearance of the 
transition evoked in the previous part. 
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Figure 9: Evolution with the vacancy density of \s, analogue of a susceptibility. 



4 Phase diagram 

4.1 Representation of the phase diagram 

We are now in position to build the full (p, T) phase diagram. We have seen that the system can end 
up in four possible states: (i) segregated state, (ii) mixed state, (iii) diluted segregation state, (iv) 
frozen state. The domains of existence of these 'phases', represented in the phase diagram, Fig. [l0| 
can be briefly described as follows: (1) For vacancy densities below 46% a dynamical transition line 
Tf(p) separates the frozen state from the segregated one. (2) Above this line the segregated state 
exists in a large domain bounded by a transition line T c (p). (2a) This line is almost at a constant 
value of the tolerance for p < 50%, and of a discontinuous nature (first-order like) up to p — 50%; 
it separates the segregated state from the mixed state. (2b) At p — 50%, the transition becomes 
continuous, and the line becomes almost parallel to the T axis; there it separates the segregated 
state from a diluted segregated one. (3) Beyond T c (p), at high values of the vacancy density one 
goes gradually (no sharp transition) from the diluted segregated state at low T to the mixed one at 
high tolerance values. We now discuss in more detail the various phases and transitions involved. 
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Figure 10: Phase diagram of the studied Schelling model. The blue crosses correspond to the 
continuous transition between the segregated and dilute segregated states. The red triangles and 
pluses are the upper and lower limits of the transition between the frozen and segregated states. 
The red lines separate the segregated state from the mixed one. Note that the tolerance T only 
takes discrete values. 



4.2 Transition frozen state / segregated state 

At low tolerance there is a transition where the system abruptly switches from a frozen to a seg- 
regated state. The analogues of the specific heat and of the susceptibility do not have a singular 
behavior in the vicinity of the transition (Fig. [8] and [9]). As we have seen, this is explained by the 
dynamical nature of this transition - in contrast to the other transitions which are thermodynamical- 
like. We locate the transition by looking at the jump of the segregation coefficient, both at fixed p, 
increasing T (Fig. [3| , and at fixed T, increasing p (Fig. 111. As the density of vacancies p increases , 



the transition occurs for lower values of the tolerance. When this density is higher than 46% this 
line of transition does not exist anymore. We remark that this frozen-segregated state transition 
line can fluctuate depending on the order in which the agents are chosen during the dynamics: in an 
interval of tolerance and for a given vacancy density, the system may end either in a frozen state or 
in a segregated one. The corresponding range of tolerance is given in the appendix, Table A. 4. 6 for 



each tested vacancy density. Let us notice that inside this frozen phase, any initial configuration, 
segregated or not, with randomly distributed vacancies is very close to a stationary state. 



4.3 Transition segregated state/mixed state. 

The segregated phase is upper bounded by a transition line where the clusters disappear and the 
two types of agents become mixed. 

- This change is abrupt for vacancy densities p strictly smaller than 26%. Indeed, the plots of 
the segregation coefficient against the tolerance for values of p in this range show a jump from ~ 1 
to - (Fig. J]! . 

- For 26% < p < 48%, the mixed state is reached after an intermediate state where the system 
has no dynamical stability: it can oscillate between several acceptable configurations, leading to 
large fluctuations in the segregation coefficient. 

- For a small range of values, 50% < p < 56%, the segregated states continously become mixed 
states. There, the line abruptly turns downward. This area of the phase diagram, with a change in 
the nature of the transition and a sudden downturn of the transition line, is more difficult to study 
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Figure 11: Segregation coefficient for different values of the tolerance versus the vacancy density. 



because of the discreteness of T values and possible finite size effects. 

For a given p the value of the tolerance at which the segregated-mixed transition occurs can be 
located from the position of the peak of the analogue of either the specific heat or of the suscep- 
tibility (Fig. [8] and [9]). In order to complete the study of this transition, we have considered the 
distribution of the segregation coefficient for several values of T around the transition at different 
vacancy densities (see appendix) . This illustrates the three ways to go from a segregated to a mixed 
state . 



4.4 Transition segregated state/diluted segregation state. 

As just mentioned, the boundary T c {p) of the segregated phase becomes weakly dependent on 
the vacancy density p when p becomes slightly larger that 50%. To locate the transition line we 
plot for several values of T the variation with p of the segregation coefficient C s (Fig. |TT]), and of 



the analogue of the susceptibility Xs (see appendix Fig. A.4|. At high vacancy densities, beyond 
the transition line, the nature of the phase is different at high and low tolerances. As we have 
seen, at high tolerance values there is a mixed state, whereas at low tolerance one finds a diluted 
segregation state that leads to very small segregation coefficients. Indeed, at high p (> 0.592746, 
the percolation threshold |18j). the high probability of percolation of vacancies prevents the forming 
of large clusters. In this high vacancy density domain, we do not find any sharp transition from the 
diluted segregated state to the mixed one, but only a gradual change as the tolerance increases. 



4.5 Comparison with the Blume-Capel phase diagram. 

It is instructive to compare the phase diagram with the one of the Blume-Capel model [TH [T3] 
evoked above (see appendix). For this model a transition line separates a ferromagnetic (segre- 
gated) phase from a domain where one goes gradually from a paramagnetic (mixed) phase to a 
phase where the vacancies predominate. This line changes as well its nature from discontinuous 
to continuous. However, it is the ferromagnetic-paramagnetic transition which is second order, 
whereas the segregated-mixed transition is first order-like (the change is abrupt). Conversely, the 
transition between the ferromagnetic and the high vacancy density phases is first order, whereas 
the corresponding transition in the Schclling model is continuous. 
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Conclusion 



We analyzed a variant of the Schelling model from a physical point of view. We have introduced a 
measure of segregation and analogues of physical quantities such that the fluctuation coefficient and 
the susceptibility where, remarkably, the analogy between the tolerance and the temperature proves 
fruitful. These quantities allowed to identify the different phases of the system and characterize the 
transitions between them (thermodynamical or dynamical like, discontinuous or continuous). The 
main results have been summarized as a phase diagram in the (p, T)-plane where p is the vacancy 
density and T the tolerance. Considering larger neighborhoods would allow to have a larger set of 
values for T and approach a continuous model. A more precise location of the phase boundaries, if 
needed, would require computationally costly simulations of larger network sizes. 
We have seen in particular that the segregated phase occupies a large domain (up to a tolerance T 
as high as 3/4), confirming Schelling's intuition on the genericity of the segregation phenomenon. 
The abrupt transition from a mixed to a segregated state could be interpreted as the tipping point 
- more precisely the rapid ethnic turnover - observed and studied by social scientists [55]. Besides, 
the diluted segregation state might be relevant for low-density suburban areas. The frozen state 
would probably be unstable in a more realistic model allowing for migratory flows of discontent 
agents to other cities. 

The tools and methods presented here could be used to study other Schclling-like models. Clearly 
future works should focus on models grounded on empirical data and where the agent decision rules 
take into account relevant socio-economic factors [3]. Yet it is known from a large body of work 
in statistical physics that one needs also to explore more widely the space of models in order to 
identify what makes a particular behavior specific or generic. As already mentioned our goal here 
was to provide generic tools for the analysis of models of socio-dynamics. The variant of Schelling's 
segregation model we have studied as a test of our approach has the advantage of being identical 
or very close to variants already studied in the literature and to allow links with known - and non 
trivial - spin models. 
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A Appendix 



A.l Numerical simulations 

All the simulations except on the Fig. [T] were performed on a 50 x 50 lattice. We tested all the 
meaningful values of the tolerance T at even values of the vacancy percentage p. With this choice 
for p, the number of vacancies and of agents of two colors are exactly equal to the integers p * L 2 
and L 2 (l — p). To take only even values of the vacancy percentage p also allows to moderate the 
computational cost. 



A. 2 Real-Space Renormalization procedure 



To identify clusters at a larger scale, we performed the following renormalization procedure. We 
divide the lattice into squares of 4 sites. On each of these little squares, we look at the bottom right 
site : 

- If this site and its neighborhood comprise a majority of blue (resp. red) agents, the 2x2 
square is replaced by a (single) blue (resp. red) agent. 

- If this site and its neighborhood consist of a majority of vacancies, the 2x2 square is replaced 
by a vacancy. 

- If there is no majority, the 4-site square is replaced by an agent of the same type as the bottom 
right site (or by a vacancy if that site is empty) . 

A. 3 Contact with spin-1 models 

A. 3.1 Energy of Schelling models 

For completeness we present here a correspondence, discussed elsewhere [21], between Schelling 
segregation models and spin-1 models. 

One can associate to each site i of the lattice a spin variable Cj, taking the value if the location is 
not occupied, and 1, resp. —1, for red, resp. blue, occupied sites. With these 'spin-1' variables the 
satisfaction condition for location i (including the case where i is vacant) can be written as: 

Ci ^c J + (2T-l) C4 2 ^ S 2 >0 (6) 

je(i) j'e(i) 

where j S (i) means j belonging to the neighborhood of site i. This suggests to define, as an 
analogue of the energy, 

E s = -Y,CiC 3 -KY,c 2 c 2 v (7) 

(id) 

where K = 2T - 1 (-1 < K < 1), and the index S stands for "Schelling". 

For the Schelling original model, as well as for other variants where only unsatisfied agents can 
move, one can show |21j that the energy E$ is indeed a Lyapunov function, that is a quantity which 
decreases with time during the dynamics, driving the system towards a fixed point. Note that the 
energy is not proportional to the global utility U — Ui, where Ui is 1 if agent i is satisfied, and 
otherwise. 

This function Es, is identical to the energy of the Blume-Emery-Griffiths model [11] under 
the constraint that the number of sites of each type (0,±1) is kept fixed. This spin-1 model, and 
the Blume Capel model [TU d3] corresponding to the particular case K = 0, have been used in 
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particular to modclizc binary mixtures and alloys in the presence of vacancies. In the standard 
versions of these models, the energy contains the additional term D q 2 (the sum being over 
all the sites) , so that the total number of vacancies is fixed only in average through the Lagrange 
multiplier D: 

ebeg = - c ^ - K ^ + D J2 c * 2 ( 8 ) 

The limit D — ► — oo corresponds to the absence of vacancies, i.e. the Ising model. Large positive D 
corresponds to high vacancy densities. The term D does not appear in the energy of the Schclling 
model, not because it corresponds to D = 0, but because the density of vacancies is fixed. The 
fact that Eg is a Lyapunov function for the Schelling model where only unsatisfied agents move, 
means that such a model is equivalent to a Blume-Emery- Griffiths model without thermal noise 
(zero temperature), and under kinetic constraints (e. g., no direct exchange between two agents of 
different colors is allowed). 

The standard order parameters for these spin-1 systems are the magnetization (l/-^0Si c i 
and the quadrupole moment, (l/-^)X)i c i> but here these quantities which are, respectively, the 
difference between the total numbers of agents of different colors, and the density of occupied sites, 
are kept fixed by construction in the Schelling model. . 



A. 3. 2 Blume-Capel model: phase diagram 

A striking similarity exists between the phase diagram in the (p, T) plane of the variant of the 
Schelling model studied here, and the one of the Blume-Capel model [HJ [T3] in the (D, T) plane - 
where T is the temperature and D the parameter fixing in average the vacancy density. 



The Blume-Capel phase diagram, computed for the Moore neighborhood, is shown on Fig A.l and 
should be compared with the one on Fig [10] The red part of the transition line corresponds to a 
first-order (discontinuous) transition, and the blue part to a second-order (continuous) transition. 
Below the transition line the system is in an ordered, ferromagnetic, state. Above the transition 
line, at low and medium values of D, one finds the unordered, paramagnetic, phase, and at large D 
and low temperature, a vacancy dominated phase. In the latter phase, the typical configurations 
are not strictly comparable to the ones of the corresponding phase in the Schelling model: in the 
spins models, the clusters are compact whereas it is not the case in the Schelling model. In the 
Blume-Capel model there is no frozen phase, because there is no constraint on the dynamics. 



The transition line in the Blume-Capel model, shown on Fig |A.1| has been built up by plotting 
for fixed temperature T (resp. fixed D depending on the area of the diagram dealt with) the 
magnetization versus D (resp. T) for three different sizes, and by looking at the value of D (resp. 
T) corresponding to the intersection of the three curves. 
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Figure A.l: Phase Diagram of the Blume Capel model with first and second neighours interactions 
(Moore neighborhood). The simulations were performed using the Heat Bath algorithm. 



A. 4 Complementary analysis of the segregation model 

A. 4.1 Density of satisfied agents 

The density of satisfied agents is a good indicator of the convergence of the system towards its 



stationary state (whenever it exists), in which the fraction of satisfied agents is constant. Fig. A. 2 
shows the evolution of the density of satisfied agents for a case with low vacancy density p = 5% 
and moderate tolerance T — 0.5. We note that, in the unstable part of the phase diagram (close to 
T c at 26% < p < 50%), the density of satisfied agents is stationary: this shows that, more generally, 
the convergence of the fraction of satisfied agents does not guarantee that the system itself has 
reached a steady state. 




Figure A. 2: Evolution of the proportion of satisfied agents versus the number of steps for p = 5% 
and T = 0.5. About 60% of the agents are initially satisfied. The dynamics quickly allows the 
agents to be almost all satisfied, after only 10 steps the density of satisfied agents is very close to 
1. It increases slowly afterwards, with small fluctuations due to the possibility for satisfied agents 
to keep moving. 
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A. 4. 2 Unwanted vacancies 



The variation with the tolerance of the unilaterally unwanted density of vacancies p (Fig |A.3[ ) con- 
firms that the mixed situation observed for low values of T is due to the dynamics. The agents 
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Figure A. 3: Density p of vacancies where no type of agent would be satisfied, for several values of 
the vacancy density. 



reject all the empty spaces, consequently the system cannot evolve. At tolerances corresponding to 
the frozen state - segregated state transition, the situation reverses. All vacancies are acceptable 
for at least one type of agents. 



A. 4. 3 Analogue of the susceptibility 

The analogue of the susceptibility for the Schelling model is given by: 

I 2\ I \2 



Xs 



T 



(9) 



where the notation ( ) means an average over the configurations taken by the system after reaching 
equibrium, and T is the tolerance. These fluctuations allow to locate the transition at fixed tolerance 




(0 
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Figure A. 4: Analogue of the susceptibility for different values of the tolerance versus the vacancy 
density. The averages have been computed on 30000 simulations after equilibrium. 
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T as the density of vacancies increases. 



A. 4. 4 Blume-Emmery-Griffiths energy vs Blume-Capel energy 

We have shown that the appropriate energy related to the Schelling model is the Blume-Emmery- 
Griffiths energy at constant number of sites of each type (red and blue agents, and vacancies), 
defined by Q. However we find that the analysis of the transitions can be done as well from the 
Blume-Capel energy, that is here, since the number of sites of a given type is fixed (no D-term): 

Ebc = ~ c ^ ( 10 ) 

(id) 

Note that in the absence of vacancies, c, = ±1 so that Ebc would reduce to the standard Ising 
energy. One finds that the fluctuation coeffi cient obtained with the Blume-Capel energy, defined by 
C' s = \BaoL^sm and shown on Fig A.5 is very similar to the one obtained from the Schelling 



energy Es- Actually, in the evolution of the mean of the total energy Es, shown on Fig. |A.6[ one 
recognizes the contribution from the Blume-Capel energy Ebc (shown on Fig. [7]), simply increased 
by an additional term linear in K = 2T — 1 . These observations can be explained by the fact that 
the difference between these two energies is proportionnal to the numbers of pairs of occupied sites 
of which the fluctuations are weak. Indeed, one observes that the vacancies remain approximatively 
uniformly distributed. This comes from the fact that each agent searches for locations where the 
number of unlike neighbors is inferior to a given proportion of the number of neighbours. Note that 
the vacancies would not be uniformly placed if each agent wanted less than a fixed number of unlike 
neighbors, as considered in [T4] , 




Figure A.5: Variation of the analogue of the 
specific heat obtained from the Blume-Capel en- 
ergy for different values of the vacancy density. 



Figure A. 6: Variation of the mean of the energy 
Es for different values of the vacancy density. 
The data are normalized by 4L 2 (1 — p). 



A. 4. 5 Segregation coefficient 



In order to get more information about the segregated state-mixed state transition, it is instructive 
to look at the distribution of the segregation coefficient < s > in the vicinity of the transition (Figs 
A. 7 and A.8l. The distributions have been obtained from 30000 calculations of the segregation 



coefficient . 

At vacancy density p lower than 26%, the transition is well marked. This distribution is centered 
near at the transition point whereas it is centered near 1 at the closest inferior value of the 
tolerance. Indeed, the distributions of the segregation coefficient for successive tolerances around 
the transition are clearly separated for the vacancy density p = 24% (left Fig. A.7). 



As p increases, the transition is achieved via an intermediate state. One distinguishes three kinds of 
distributions (Figs. A.7 right and |A.8| left) which are centered around a very small value (~ 0.2), 
peaked close to 1, or centered around an intermediate value. The latter case corresponds to a 
broader distribution. 
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Figure A. 7: Distributions (normalized by the number of measures) of the segregation coefficient for 
p = 24% and p = 28%. 
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Figure A. 8: Distribution of the segregation coefficient for p = 44% and p — 50%. 



Once the vacancy density is greater than 50%, all the distributions of the segregation coefficient 
begin to blend together (see Fig A. 8 right). This confirms that the increase in the vacancy density 
is accompanied by a change in the nature of the transition between 46% and 50%, which becomes 
continuous. 
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A. 4. 6 "Frozen-segregated" transition line 



The transition line between the frozen and the segregated states has been determined by locating 
the jump of the segregation coefficient < s > from ~ to ~ 1. The initial configuration and the 
order of choice of the agents, when the dynamics is applied, create fluctuations on the limit between 
the two states. Actually, for some sets of parameters (p, T), the system may end either in blocked 
or in segregated configurations depending on the order of choices of the agents. One cannot exclude 
that this unstable domain is due to finite size effects, hence disappearing in the infinite network 
size limit. However, it is not surprising to find such a metastability effect close to a discontinuous 
transition. 
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Table 1: "Frozen state-segregated state" transition line. The limits between which the system may 
end in a frozen state or segregated state depending on the order of the dynamics are obtained by 
performing 100 tests on which we look at the percentage of "frozen states". If for given T and p, 
the percentage of frozen states (resp. segregated states) is very high (> 95%), we consider that the 
corresponding equilibrium configuration is a frozen one (resp. segregated). 
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